We investigate the moments of a smooth counting function of the zeros near the central point of L-functions of weight k cuspidal newforms of prime level N. We split by the sign of the functional equations and show that for test functions whose Fourier transform is supported in (-1/n, 1/n), as N --> oo the first n centered moments are Gaussian. By extending the support to (-1/n-1, 1/n-1), we see non-Gaussian behavior; in particular the odd centered moments are non-zero for such test functions. If we do not split by sign, we obtain Gaussian behavior for support in (-2/n, 2/n) if 2k >= n. The nth centered moments agree with Random Matrix Theory in this extended range, providing additional support for the Katz-Sarnak conjectures. The proof requires calculating multidimensional integrals of the non-diagonal terms in the Bessel-Kloosterman expansion of the Petersson formula. We convert these multidimensional integrals to one-dimensional integrals already considered in the work of Iwaniec-Luo-Sarnak, and derive a new and more tractable expression for the nth centered moments for such test functions. This new formula facilitates comparisons between number theory and random matrix theory for test functions supported in (-1/n-1, 1/n-1) by simplifying the combinatorial arguments. As an application we obtain bounds for the percentage of such cusp forms with a given order of vanishing at the central point.